Kinesin as an electrostatic machine 
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Kinesin and related motor proteins utilize ATP fuel to propel themselves along the external surface 
of microtubules in a processive and directional fashion. We show that the observed step-like motion 
is possible through time varying charge distributions furnished by the ATP hydrolysis circle while 
the static charge configuration on the microtuble provides the guide for motion. Thus, while the 
chemical hydrolysis energy induces appropriate local conformational changes, the motor translational 
energy is fundamentally electrostatic. Numerical simulations of the mechanical equations of motion 
show that processivity and directionality are direct consequences of the ATP-dependent electrostatic 
interaction between the different charge distributions of kinesin and microtubule. 



I. INTRODUCTION 

The kinesin family is a set of motor proteins that move 
on the surface of microtubules and shuttle various cargo 
molecules to different parts of the cellQ-0. The energy 
for kinesin transport motion is provided by hydrolysis of 
ATP molecules that are in its vicinity and attach to a 
specific site of the protein. Despite much progress in the 
area, the precise, microscopic fashion of ATP action is 
not known. Kinesin motion is step-like with one step for 
each ATP hydrolysis while processivity and directional- 
ity depend on specific neck properties of each particu- 
lar type of the family. These features of motor protein 
motion led several years ago in the introduction of me- 
chanical ratchet models that could provide some insight 
on the phenomenology of individual 0, 0, 0, Q, [u| 
and collective 0, 01 motion. Although many issues 
and especially those related to the random aspects of 
the walk where addressed successfully by these models, 
fundamental questions on the energetics and nature of 
forces that enable the walk have not been discussed. In 
the present work we focus on the latter; specifically we 
use electrostatic information of the microtubule and ki- 
nesin and evaluate the full electric force that binds the 
latter on the former. We then show that the actual ki- 
nesin walk is caused by charge distribution changes on 
kinesin enabled by the action of ATP. The walk that is a 
combination of conformational changes accompanied by 
kinesin charge reshuffling is found to be fundamentally 
governed by electrostatic forces. 

II. KINESIN ELECTROSTATIC MODEL 

Kinesins walk on microtubules; the latter are made 
by alpha and beta of tubulin proteins that are highly 
charged. Molecular modeling shows that, in the presence 
of water, tubulins have additionally a permanent dipolc 
moment that points towards the cylindrical symmetry 
axis of the microtubule with a non-zero component par- 
allel to the protofilaments towards their minus endQ. 



Electrostatically thus a microtubule has negative surface 
charge accompanied by a positive charge distribution 
in its immediate interior with a slanted polarization 
vector that is larger in the beta than the alpha subunit 
(Figure 0. This symmetry breaking induced by the 
tubulin dipolc moment direction plays a key role in the 
specificity of the attachment of a kinesin head onto a 
beta tubulin subunit 2]. 

Structural analysis of kinesin and related proteins sug- 
gests that there are three charge domains, one in each 
head region and a third one in the central neck linker do- 
main. The head charge distribution depends strongly on 
the ATP or ADP presence while the neck charge is spe- 
cific to the protein. For modeling purposes it is sufficient 
to consider kinesin and similar motor proteins as a body 
containing three charges located on the apexes of a trian- 
gle with angle 7 subtending from the neck charge to the 
two head charges (Figure ^Jd) These charge distributions 
interact with tubulin multipole fields and determine at 
each instant of time the dynamic state of the molecule. 

In our analysis of kinesin walk to be detailed below, 
we place motor proteins and microtubules in an over- 
damped medium with thermal energy ksT = A.lpNnm 
and motor drag force 10~ 6 pNs/rad per rotational de- 
gree of freedom. For the sake of simplicity, we assume a 
relative electric permitivity equal to 80 everywhere and 
a Debye length Id above 3nm. Although Dcbye-Hiickel 
theory predicts Id ~ 1, the interaction range may be en- 
hanced by the channeling of the electric field along the 
interior of the proteins, which is not accesible by diffusing 
ions or water molecules. Then it is reasonable to take Id 
of the order of the size of the kinesin. We take consec- 
utive tubulin unit distance to be 8nm, Anm per a or (5 
subunit. The lateral space between two protofilaments is 
5nm while the pitch of lnm between consecutive protofil- 
amens is also included. For a processive plus-ended ki- 
nesin and its non-processive-minus-ended chimera ncd we 
use for the head-to-hcad distance the value of 6nm while 
we consider the heads to be charged with +2e in the ab- 
sence of nucleotides, — 2e with ATP attachment and — le 
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FIG. 1: a) Microtubule electrostatic model with arrows indi- 
cating local dipole moments. The a tubulin subunits (dark) 
have smaller dipole moment that the f3 units (lighter color), 
b) Electrostatic configuration of kinesin and tubulin made 
protofilament prior to ATP hydrolysis. Angle 8 is polar while 
4> is azimuthal. The central-neck charge sign depends on the 
type of molecular protein while the head charges depend on 
the ATP hydrolysis circle. Dipolar lenghts d a , dp and dipolar 
angles ui a , uip are different in a and (3 subunits respectively. 



with ADP. These head charge values are consistent with 
the charges involved in hydrolysis reaction 



ATP 



OH 1 



ADP* 



H 2 POi 



(1) 



The value of the central motor charge changes for dif- 
ferent proteins; as a result we consider different motors 
with charges in the range [— 2e : +2e]. 

For the electrostatic distribution of microtubules we 
assign negative surface charge q = —27 e per tubu- 
lin subunit while include a positive charge distribution 
in the interior leading to dipole moment magnitude of 
p = 5000Z? ~ 100Cnm[2, or d ~ 4nm (p = qd). Finally, 
we use a dipolar tilt with length and angle values equal to 
d a ~ 2nm, u> a ~ 0.07rad and dp ~ 4nm, up ~ O.lArad, 
for the a and (3 subunits respectively (Figure H 3 ). 

The microtubule-induced kinesin interaction potential 
is given by 



N 
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where fi is the position vector labeling the charges on 
kinesin (i — 1,2,3), while rj is the location of the N 



microtubule charges qj on the a and /? subunits. k is 
the inverse of the Debye length, which we take around 
~ 3.5nm, and a ~ Inm is the excluding volume radius 
as described in Ref. |18| . A regime with a greater a and 
lower Debye length Ijj is also operative. We considered a 
flat microtubule with five protofilaments; due to the rapid 
decay of the force out from protein volumes, we include in 
total the N — 10 closest tubulin charges to kinesin. For 
the simplest case when neck and head charges are aligned 
we have 7 = n and the protein reduces to a triply charged 
rigid rod. 

We have tested the model for several values of 7 smaller 
than 7r; we have found that the model is not fully opper- 
ational for 7 < 175°. We have also tested the model 
considering that the angle 7 has some elasticity as well. 
We found that in case the angle is very soft, the raising 
head rapidly loses its power to pull the motor forward. 
For the purposes of the model, the motor core and the 
neck need to be stiff enough to maintain the values of 7 
in the range discussed above even at maximum load con- 
ditions. In this simpler rod-like configuration the polar 
angle 9 and the azimuthal angle 4> are sufficient for de- 
scribing the motor rotation. We can describe the motion 
using the following overdamped equations of motion for 
kinesin: 
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where A is the drag coefficient, L the head-to-head dis- 
tance and V(r) the total microtubule electrostatic po- 
tential of Eq. J5J at the Cartesian location r. The en- 
vironment is simulated through the thermal forces £e(t) 
and £0; for each we have (£(*)) = and = 
2\k B T6(t — t'). In order to integrate Eqns. (21 0J we 
need to perform at each instant of time the Cartesian- 
to-polar transformation 6 = 9(f, L). For the rotation we 
consider the attached head to coincide with the origin 
of the coordinate system; the latter is shifted by 8nm 
each time a step is completed. The simple Larmor-like 
rotation of the protein for 7 = tt becomes a more com- 
plex rigid body rotation for 7 < tt. Extensive numerical 
simulations with the model described previously for the 
kinesin-microtubule complex led to the following quan- 
titative picture for the motion. In equilibrium (parked) 
state the nucleotide-free head is positively charged and 
thus attached to the beta subunit of the negative micro- 
tubule surface while the other head containing ADP is 
negatively charged and tethered through the neck. Due 
to the slanted microtubule dipole moment, the kinesin 
axis is tilted but in a direction determined by the cen- 
tral charge sign; for wild type kinesin points to the fast- 
growing-end while for ncd points to the minus end ( Fig- 
ure El- 



3 



a) b) 




FIG. 2: Electrostatics driven motor walk: Wild-type kinesin 
(+ states) and ncd (- states) stepping process (left column) 
and numerically determined binding protein-microtubule elec- 
trostatic potential as a function of local polar angle 9 (right 
column). (0) Before ATP hydrolysis, both kinesin (positive 
neck) and ncd motors (negative neck) are in parked config- 
uration pointing in opposite directions due to the difference 
in the central charge. The corresponding equilibrium angles 
are determined for the minimum potentials. (+1) ATP entry 
in the kinesin (attached) head pocket, with an accompanying 
charge change while the ADP at the other (tethered) head be- 
comes unstable. (+2) The reversal and shift of the interaction 
potential of the previous state leads to falling of the tethered 
head deterministically towards the plus end. Since the length 
of the motor is not sufficient for reaching the next tubulin 
subunit we have (+3) a detachment and rising of the trail- 
ing head in such a way that allows the other head to slide to 
the next binding site. The ncd motor protein cycle proceeds 
similarly (- states) but the parked state is tilted towards the 
minus-end. Moreover, the negative-central charge induces a 
potential shift which is opposite to the plus-ended case. The 
falling of ncd motors is slower than positive-charged-neck mo- 
tors, so the probability that the attached heads begins the ris- 
ing before the tethered heads completes the falling is greater, 
leading to non-processivity. 



All experimental data agree that after nucleotide en- 
try in the head pocket, the binding of the attached head 
to the microtubule and the ADP binding to the tethered 
head become unstable, signaling the onset of a cascade 
of events leading to the walk. In our model this criti- 
cal juncture occurs because in the absence of nucleotide 
the positively charged head becomes negative acquiring 



— 4e charges after ATP binding. Although the entry of 
a negatively charged molecule in the head pocket is nec- 
essary for the commencement of the walk process, the 
release of the attached head is not immediate, since a 
head-tubulin chemical bond needs to be broken. ATP 
hydrolysis energy thus is used for a local conformational 
change that captures the negative charge after dephos- 
phorylation. Although charged ADP nucleotides in the 
medium compete with ATP entering also in the pocket, 
they do not hydrolyze and, as a result, cannot induce 
the local conformation that will trap them. In our elec- 
trostatic model thus hydrolysis energy induces primarily 
a specific local conformational change, being the power 
stroke electrostatic in nature. 



III. KINESIN DYNAMICS 

Based on this picture, the protein parked configura- 
tion is a result of the microtubule-empty head attrac- 
tion, microtubule-tcthered head repulsion while the ki- 
nesin angle wrt microtubule is due to the tubulin dipolar 
tilts. This stable configuration is maintained during a 
random and [ATP]-load-dependent dwell time until a new 
ATP binds into the attached head. When this occurs, the 
local charge changes and the ADP trapped in the teth- 
ered head experiences an additional repulsion, becoming 
unstable and eventually opening the pocket and exiting. 
After ADP expulsion the tethered head becomes positive, 
the attached head charge is negative and the protein be- 
comes electrostatically unstable. As a result the tethered 
head collapses onto the microtubule ("falling processes"), 
while the attached head is forced to detach from the sur- 
face ("rising process") and a new, shifted, park state is 
reached. 

Every kinesin step thus includes falling of the lead- 
ing head onto the microtubule followed by rising of the 
trailing head ( Figure [5J . We performed extensive 3D 
numerical simulations and found that the electrostatic 
force field generated by the microtubule is able to drive 
these two sequential processes while keeping the motor 
faithful to a given protofilament. We note that, energeti- 
cally, the parked state represents a potential energy min- 
imum for the interaction between free head charge and 
microtubule. When ADP is released and the electrostatic 
charge distribution of the motor changes this minimum 
becomes a potential maximum, i.e. the repulsion of the 
tethered head turns into attraction. 

For a null-charged-neck case, even though the motor is 
tilted in the parked state, there is no motion direction- 
ality since the falling of the tethered head does not have 
a preferred collapsing side. The directionality features of 
the motor enter through the charge distributions of the 
neck region, which is responsible of the potential shift 
between stable and unstable equilibrium points. Specif- 
ically, when the central charge is positive, not only the 
protein is more processive due to the attraction with the 
negatively charged tubulin surfaces, but it also walks to- 
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wards the plus-end of the microtubule; this is precisely 
the case with wild-type kinesin. In ncd, on the other 
hand, the central charge is negative, the protein is non- 
processive and walks in the opposite direction towards 
the minus end. Additionally, in proteins with no charge 
in the neck region, the walk is known not to be deter- 
ministic but random 3] . We note further that our simu- 
lations are compatible with the experiments in reference 
whereby the falling motor direction is parallel to the 
protofilament axis, with some variations and latteral col- 
lapses leading to protofilament changes due to thermal 
fluctuations. The lateral periodicity of the microtubule 
lattice is 5nm while the axial one only 4nm; these values 
determine crucially the motion after the destabilization 
of the parked state since the electric force is stronger in 
the axial direction. We find that this electrostatic model 
captures fully the directionality and processivity features 
of the known motors. 

Two different time-scales characterize the actual ki- 
nesin motion. While dwell times between two steps, mod- 
ulated by ATP concentration and the external load are 
of the order of milliseconds Q, the action of the step itself 
lasts only some microseconds Q ; this is the reason for the 
stepping appearance of the trajectories. In the relatively 
short period of the step time, several processes must oc- 
cur. First, falling of the leading head after ADP release, a 
process activated by ATP binding on the attached head. 
Subsequently the two heads remain attached onto the mi- 
crotubule for a time not longer than 30 microseconds . 
Recent measurements were able to track bead-movements 
at the microsecond scale 6]. The recording of a single 
step, in our interpretation, reveals two different regimes 
that may be separated by a quick stage at which both ki- 
nesin heads are attached to the microtubule. Kinesin cry- 
tallization data show that the separation between heads 
is smaller than 8 nmQ, i.e. slightly smaller than the 
protofilament spatial period and, as a result, after falling 
the leading head cannot reach the beta subunit location 
before the trailing head detaches. Thus, sliding may take 
place while the trailing head rises to reach the new parked 
state. This aspect is portrayed in Figure|3]in a sequential 
form although in practice both sliding and rising could 
occur simultaneously. 

In order to address motor protein processivity we 
consider the two sets of chemical reactions that arc 
activated by ATP binding at the trailing head. The first 
group is related to a cascade of rapid reactions involving 
ATP hydrolysis, dephosphorylation and subsequent head 
detachment with an estimated time 500 /is for eachQ. 
Secondly, ATP activates ADP release in the other head, 
although the typical time of this process is not known. 
There is a competition between these two reactions; 
if hydrolysis and detachment occur faster, the whole 
motor detaches from the microtubule and processivity 
ceases. If, on the other hand, ADP release and falling 
is the fastest process, then a small substep occurs while 
both heads are attached. Processivity thus is related 
to the competition of these two processes. In the case 
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FIG. 3: Numerical results for a kinesin step at (a) ms time- 
scale and (b) fis time-scale, for 7 = ir; s-axis denotes time 
while i/-axis kinesin position.. Occurrence of non-processivity 
and substeps are in competition since for the former detach- 
ment of the trailing head occurs before than the falling of 
tethered head, while for substeps to occur the falling process 
must be done before detachment, leading to a processive mo- 
tion. The millisecond time-scale numerical results shown in 
(a) are detailed at a microsecond time-scale in (b) where the 
aforementioned competition is seen. 



of non-processive ncd, for instance, the existence of a 
negative central charge has the effect of slowing down 
the falling process while speeding up the rising one 
resulting, thus, in a non-processive motor. This effect 
has been tested experimentally in Ref. [l9|. 



We point out that when ATP enters into the at- 
tached head, the system becomes electrostatically un- 
stable. However, detachment of the head occurs some 
time later triggered by the phosphate release. This fact 
implies that the chemical bonding that binds the motor 
to the microtubule, while stable in absence of ATP, be- 
comes unstable when the nucleotide arrives. The time 
delay between ATP entry and bond dissociation appears 
to be crucial for motor processivity. In the Figure [21 we 
have assumed that both possesive and non-possesive mo- 
tors may perform the full detachment-attachment cycle, 
although, in the case of ncd, the chemical bond may be 
actually broken before the falling head reaches the next 
binding site. The time duration of a single kinesin step is 
very small compared to the parked dwell times, and, as 
a result, global observables such as the mean velocity or 
randomness can be predicted using only chemical kinet- 
ics. This has been shown in reference 20] where a kinetic 
approximation was shown to be sufficient in fitting the 
measured data, even in the presence of an external load. 
However, our present approach shows, that thermal fluc- 
tuations alone may not produce directly the 16nm long 
displacements in a jis time-scale necessary for the walk. 
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IV. CONCLUSIONS 

In this work we attempted to analyze the complex mo- 
tion of molecular motor proteins from a mesoscopic point 
of view placing emphasis on the fundamental interac- 
tions that enable the motion. We found that the mo- 
tion is primarily driven by the electrostatic interaction 
between the charged microtubule surface and the fluctu- 
ating motor head charges. The nonequilibrium aspect of 
the walk is provided by the ATP hydrolysis cycle that fur- 
nishes appropriate charge motor distributions that make 
the walk possible. We made several assumptions in this 
work, most of which have been motivated by the current 
status of knowledge in the area. We considered the mo- 
tor as a relatively rigid one; this assumption is presently 
well founded, however it may be lifted through an im- 
proved motor model that includes additionally protein 
eleasticity. Furtheremore, our model docs not address at 
all the structual changes in the head pockets effected by 
ATP hydrolysis, the mode of the local energy relase as 
well as the mechanism for the ADP detachment from the 
thethered head. These issues that may involve complex 
conformational changes as well as possibly charge or en- 
ergy transfer processes must be addressed through a more 
foundamental, microscopic model. Finally, it is crucial to 
evaluate the range of Debye screening in the vicinity of 
a highly charged surface such as the microtubule. While 
our motor model remains fully operational to a Debye 



length range of approximately 3nm, one certainly needs 
to address the complex nature of the electrostatic shield- 
ing in the vicinity of the microtubule and assess the true 
range of the electrostatic forces involved in the walk. We 
note that our model provides a consistent dynamical pic- 
ture for the kinesin walk that is based on two premises, 
viz. (a) that the ATP hydrolysis energy is used for a 
head local conformational change that captures locally 
charges and (b) the motor motion is driven by electro- 
static forces. The emerging qualitative and quantitative 
picture for the walk is fully compatible with all known 
experimental data, while, furthermore, it is testable ex- 
perimentally. If electrostatics, charge transfer as well as 
capture indeed link the chemistry of ATP and the me- 
chanics of kinesin, this may hold true for other ATP- 
dependent processes as well. 
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